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Lattice QCD with background magnetic fields is used to calculate the magnetic moments 
and magnetic polarizabilities of the nucleons and of light nuclei with A < 4, along with the 
cross-section for the Ml transition np —> d' y, at the flavor SU(3)-symmetric point where the 
pion mass is m n ~ 806 MeV. These magnetic properties are extracted from nucleon and 
nuclear energies in six uniform magnetic fields of varying strengths. The magnetic moments 
are presented in a recent Letter [T] . For the charged states, the extraction of the polarizability 
requires careful treatment of Landau levels, which enter non-trivially in the method that is 
employed. The nucleon polarizabilities are found to be of similar magnitude to their physical 
values, with = 5.22 ( ) (0.23) x 1CT 4 fm 3 and /3 n = 1.253 ( 1™®® ) (0.055) x 1CT 4 fm 3 , 

exhibiting a significant isovector component. The dineutron is bound at these heavy quark 
masses and its magnetic polarizability, (3 nn = 1.872 ^ ^ (0.082) x 10 -4 fm 3 differs 

significantly from twice that of the neutron. A linear combination of deuteron scalar and 
tensor polarizabilities is determined by the energies of the j z = ±1 deuteron states, and 
is found to be Pd.,±i = 4.4 ( ® ^ (0.2) x 10 -4 fm 3 . The magnetic polarizabilities of the 

three-nucleon and four-nucleon systems are found to be positive and similar in size to those 
of the proton, /33 h 6 = 5.4 ^ 43 ’i ) (0-2) x 10 -4 fm 3 , fan = 2.6(1.7)(0.1) x 10 -4 fm 3 , /34 He = 

3-4 ( 1^9 ^ (0.2) x 10~ 4 fm 3 . Mixing between the j z = 0 deuteron state and the spin-singlet 
np state induced by the background magnetic field is used to extract the short-distance two- 
nucleon counterterm, Lx, of the pionless effective theory for NN systems (equivalent to the 
meson-exchange current contribution in nuclear potential models), that dictates the cross- 
section for the np —> dy process near threshold. Combined with previous determinations 
of NN scattering parameters, this enables an ab initio determination of the threshold cross- 
section at these unpliysical masses. 


PACS numbers: 11.15.Ha, 12.38.Gc, 13.40.Gp 
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I. INTRODUCTION 

The charge, magnetic moment and magnetic polarizability of a composite system describe its 
linear and quadratic response to a uniform, time-independent magnetic field. These properties are 
determined by the distribution of the constituents of the system and by the currents induced by 
the field. As such, measurements of the magnetic properties of the nucleons and nuclei provide 
important information about their internal structure. Furthermore, these quantities also serve 
to parameterize the cross-section for low-energy Compton scattering on such targets. Histori¬ 
cally, the magnetic moments of nucleons and light nuclei provided some of the first indications 
of substructure, and by now they are well known. The primary focus of this article is on the 
magnetic polarizabilities and the cross-section for the radiative capture process np —y d'y at low 
energies which is dominated by the Ml multipole. While the magnetic polarizability of the pro¬ 
ton 1 , f3p Xp = (3.15 ± 0.35 ± 0.20 ± 0.30) X 10 -4 fm 3 , is well determined experimentally J3HS] , the 
magnetic polarizability of the neutron /3^ p = (3.65 ± 1.25 ± 0.20 ± 0.80) x 10 -4 fm 3 remains quite 
uncertain 130 M]- This uncertainty is largely a consequence of the lack of a free neutron target; 
the neutron polarizability must be determined from that of light nuclei, primarily the deuteron 
(see Ref. m for a recent summary). The smallness of the nucleon polarizabilities, compared with 
their spatial volumes, ~ 1 fm 3 , indicates that they are quite magnetically rigid, with the spins and 
currents of their constituents influenced only at a modest level by external fields. The uncertainty 
in /3f xp is large enough so that a significant isovector polarizability remains a possibility. For a 
recent review, see Ref. [12]. 

From a theoretical standpoint, the leading contributions to the nucleon magnetic polarizabilities 
result from both pion-loop effects and A-resonance pole contributions J32H32]- The A-pole contri¬ 
bution is 0(e 2 /[M^(M/\ — Mtv)]) and is considerably larger than the experimentally determined 
polarizabilities. At the physical light-quark masses, significant cancellations between the param¬ 
agnetic (A-pole) and diamagnetic (loop) contributions are in effect. As the various contributions 
depend differently on the light-quark masses (the leading pion-loop contributions scales as ~ m~ l 
while the A-pole contribution is only weakly dependent on the quark masses), it is expected that 
(dp and f3 n will vary reasonably rapidly near the physical point. 2 Because it is relatively mass 
independent, the A-pole contribution provides an estimate of the expected size for polarizabilities 
at any quark mass and it will be used to assess the naturalness of the polarizabilities extracted 
from the present Lattice QCD (LQCD) calculations. 

The leading contributions to the magnetic polarizabilities of light nuclei are determined by the 
nucleon charges, magnetic moments and polarizabilities. At a sub-leading level, the forces between 
nucleons, both those responsible for nuclear binding and meson-exchange currents (MECs), modify 
these one-body contributions and produce short-distance contributions that are unrelated to the 
electromagnetic interactions of single nucleons. Accounting for this in a consistent manner is 
non-trivial and requires a well-controlled power-counting. This has been investigated in detail 
for the deuteron m and light nuclei [20]. Experimentally, the polarizabilities of the deuteron 
have been measured through Compton scattering (leading to extractions of the neutron magnetic 
polarizability) and further measurements will be performed with increased precision at the HlyS 
facility [21], MAX-Lab at Lund [22] and at MAMI in Mainz |23j . Plans for a new generation of 
Compton scattering experiments on other light nuclei are also being developed (see, for example, 
Refs. HH). 

The radiative capture process, np —y dy, and the inverse processes of deuteron electro- and 


1 The experimental polarizability reported here is defined with Born terms subtracted; the total 0(B 2 ) shift in the 
energy of the proton is larger by (167rM 3 ) -1 = 0.15 x 10 _4 fm 3 [2]. 

2 Recent chiral effective field theory calculations support this expectation 118] . 
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photo-disintegration, 7 —> np , are important in early universe cosmology and have led to crit¬ 
ical insights into the interactions between nucleons, and in particular, interactions with photons. 
At very low energies, the Ml magnetic multipole transition amplitude is dominant and is de¬ 
termined primarily by the isovector nucleon magnetic moment. The short distance contributions 
(equivalently, MECs) are sub-leading and modify the cross-sections at 0(10%) [ZD, 28j and are well 
determined by experiment in this particular case. Given the phenomenological importance of this 
and other related processes, it is important to understand these contributions from first principles. 

In this work, we present lattice QCD (LQCD) calculations of the magnetic moments and po¬ 
larizabilities of the proton, neutron and s-shell nuclei up to atomic number A = 4. Further, the 
jz = I z = 0 np systems are used to investigate the short-distance two-nucleon contributions to 
the cross-section for np —> d'y. The methods and calculations that are presented are an extension 
of those used to calculate the magnetic moments of light nuclei in Ref. |T] and the result for the 
capture cross-section have been recently highlighted in Ref. [ 25] . In Section [TTJ the methodology of 
using background magnetic fields in LQCD calculations to determine hadron magnetic moments 
and polarizabilities and magnetic transition amplitudes is presented. Section III discusses the 
results of our calculations, by first discussing the general analysis methods used to extract the 
magnetic properties from energy shifts before turning to a discussion of each of the nucleons and 
nuclei that are studied. Our results are summarized in Section [TV] , and opportunities for possible 
extensions of this work are also outlined. Finally, Appendix [A] is dedicated to defining the for¬ 
malism underpinning the LQCD methodology used in the present calculations. Explicit examples 
centered on the behavior of point-like charged particles are explored. 


II. METHODOLOGY 

In order to calculate the magnetic structure of nucleons and light nuclei, and the low-energy 
cross-section for np —> dj, it is sufficient to perform LQCD calculations of these systems in uniform 
time-independent magnetic fields. In sufficiently weak background fields, the quantities of interest 
can be extracted directly from the energy eigenvalues of the spin states of these systems. 


A. Background Magnetic Fields 


In this work, lattice QCD calculations are performed using one ensemble of QCD gauge-field 
configurations with Nf = 3 degenerate dynamical flavors of light quarks on a L 3 x T = (32 3 x 48) a 4 
discretized spacetime. This ensemble was generated using a Liischer-Weisz gauge action m and 
a tadpole-improved clover-fermion action m at a coupling corresponding to a lattice spacing of 
a = 0 . 110 ( 1 ) fin |321 (further details of this ensemble can be found in Refs. [331 HU)- All three 
light-quark masses in this ensemble are equal to that of the physical strange quark, producing a 
pion of mass ~ 806 MeV. In the present set of calculations, OQO 3 ) gauge-field configurations 
separated by 10 Hybrid Monte-Carlo trajectories were used from this ensemble. 

Background electromagnetic fields have been used extensively in previous calculations of the 
electromagnetic properties of single hadrons, such as the magnetic moments of the lowest-lying 
baryons mmm, the polarizabilities of mesons [381 HD \MM and the electric polarizabilities of 
baryons |46L 147] . In addition, the magnetic moments of light nuclei have been recently calculated 
in Ref. [lj. In order for the quark fields, with electric charges q q ( q u = +| and qd )S = — | for the 
up-, down- and strange-quarks, respectively), to encounter uniform magnetic flux throughout the 
lattice, the field strength is quantized according to the condition [48] 


eB = j 2 ne z 


n € Z, 


( 1 ) 
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where e is the magnitude of the electron charge and we consider explicitly fields in the X 3 = z 
direction as indicated by the unit vector e z . Neglecting the electric charges of the sea-quarks, a 
quantized, time-independent and uniform background magnetic field oriented in the positive X3- 
direction can be implemented by multiplying the QCD gauge link variables by Uq(\) link fields, 
of the form 


u[ Q \x) 

U?\x) 

U?\x) 

Uf\x) 


exp [-iq q h^^ 


exp 


iq q n- 


L 

6ira x\ 
~LP~ 


for x\ ^ L — a 
for xi = L — a ’ 


1 , 
1 , 


( 2 ) 


where h is an integer satisfying |h| < ^L 2 /a 2 (in Euclidean space, x 4 corresponds to the Wick- 
rotated time coordinate). In physical units, the background magnetic fields used in this work are 
quantized in units of eB ~ 0.059 h GeV 2 e z ; in comparison to the nucleon mass, the dimensionless 
ratio e|B|/Mjy ~ 0.013 for the smallest field suggesting the deformations arising from the magnetic 
field are perturbatively small compared to QCD effects for \n\ < 10. As m u = rn^ = m s in these 
calculations, the up-quark propagator in the h = 1 field is the same as the down- or strange-quark 
propagator in the h = —2 field, with similar relations for the other magnetic field strengths. To 
optimize the re-use of light-quark propagators in the calculation, different quark electric charges 
were implemented by using a different magnetic field (with the same charge). Consequently, the 
Uq(1) fields that are used in this work correspond to n = 0, +1, —2, +3, +4, —6, +12, corresponding 
to magnetic fields of eB ~ 0, 0.06, —0.12, 0.18,0.24, —0.36,0.71 GeV 2 e z , respectively. 

In the presence of a time-independent and uniform magnetic field, the energy eigenstates of 
a structureless charged particle with definite angular momentum along the field direction are de¬ 
scribed by Landau levels and plane waves, rather than by three-momentum plane waves alone. 
One of the subtle finite-volume (FV) effects introduced into this calculation is the loss of transla¬ 
tional invariance in the interaction of charged particles with the background gauge field. We give 
a brief description of this problem, and relegate the more technical aspects of the discussion to 
Appendix [Aj For the implementation of the magnetic held using the links in Eq. ([ 2 ]) , the lack of 
translational invariance is made obvious by the Wilson loops, 


W ] (x 2 ) 
W 2 (x 1 ) 


L/a—1 

If u[®\x + jax 1 ) = exp 
j =0 

L/a—1 

If U^\x + jax 2 ) = exp 
3 =0 


-iq q n 


, ()TTX2 


iq q n- 


L 

67TX1 

IT 


(3) 


which wrap around the x\ and x 2 directions of the lattice geometry, respectively. These exhibit 
explicit spatial dependence. Further, there are additional effects for charged-particle correlation 
functions arising from their gauge dependence. 

Because of the lack of translational invariance, the links employed in Eq. ([ 2 ]) define a spatial 
origin x = 0, where the gauge potential vanishes, A(x) = 0. In performing the present calculations, 
the source points for the correlation functions are not restricted to x = 0 but instead are randomly 
distributed within the lattice volume, approximately restores translational invariance. In the case 
of charged-particle correlation functions, this averaging leads to non-trivial effects, because the 
overlap of a given hadronic operator onto the various Landau levels depends on the source location. 
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This is explicated in Appendix [A] where additional methods of restoring translational invariance 
are discussed in the context of a structureless charged particle. 

Post-multiplication of Uq( 1) links onto the QCD gauge links omits the effects of the external 
magnetic field on the gluonic degrees of freedom through the fermion determinant. The present 
calculations therefore correspond to a partially-quenched theory in which the sea quark charges are 
set to zero while the valence quark charges assume their physical values usmsHsn. F ° r a su(3) 
symmetric choice of quark masses, as used herein, this does not affect the magnetic moments or 
the np —> d'y transition (linear responses to the field) because the Nf = 3 charge matrix is traceless 
HI 09] and couplings to sea quarks explicitly cancel at this order (indeed the isovector nature 
of the np —> d'y transition make it insensitive to disconnected contributions even away from the 
SU(3) symmetric point). However, the magnetic polarizabilities receive contributions from terms in 
which the two photons associated with the magnetic field interact with zero, one or two sea quark 
loops. The terms involving zero and two sea quark loops are correctly implemented, however the 
square of the light-quark electric-charge matrix is not traceless and terms involving the two photons 
interacting with one sea quark loop will contribute to isoscalar quantities for any non-zero charge 
matrix. In the present work, these terms are omitted for computational expediency. 3 Generally, it 
has been found that the related disconnected contributions to analogous single-hadron observables 
are small for the vector current [521 [53], and we expect that this behavior persists in nuclei. It 
is important to remember that this systematic ambiguity is restricted to the case of the isoscalar 
polarizabilities, and that the isovector and isotensor combinations, such as f3 p — /3 n , will be correctly 
calculated for the SU(3) symmetric case. 


B. Interpolating Operators and Contractions 

In order to probe the energy eigenstates of the systems under consideration, interpolating oper¬ 
ators are constructed with the desired quantum numbers. In principle, any choice of operator that 
has a overlap onto a given eigenstate is acceptable. However, poor choices will have small overlaps 
onto the state of interest and hence will give rise to “noisy” signals with significant contamination 
from other states with the same quantum numbers. For a vanishing background magnetic field, the 
energy eigenstates are also momentum eigenstates, and in order to access the ground state energy, 
it is simplest to choose interpolating operators that project onto states with zero three-momentum. 
In this work, this approach is followed for both the neutron and proton (despite the proton carrying 
a positive electric charge). The proton correlation functions are of the form 

C( P ’ S >(t; Xt ) = (0ldj, P ’ S Ht;0)O ( p S) (x u 0)l0) B , (4) 

with interpolating operators that are given by 

0 { p s \x,t) = e vk [u\( X ,t)C'y 5 d J + ( X ,t)}u k + (x 1 t), 

Of\t- p) = ^e l PV^[uV(x,t)C 75 4(x,t)]4(x,t), 

X 

6 { p S \t-,p) = ^e l P x e^[uV(x,t)C 7 5<(x,t)]4(x,t) = ^e^O^x,!), (5) 

X X 

where (.. .)b indicates ensemble averaging with respect to QCD in the presence of the Uq{ 1) links 
corresponding to a uniform background magnetic field B = Be z , and where the spin indices of the 


3 Several approaches to these terms have been considered recently mum and may be investigated in future studies 
of nuclei, although significant computational resources are required. 
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operators, carried by the third quark, are suppressed. In Eq. ([5]), q(x) corresponds to a quark field of 
flavor q = u,d that has been smeared (54) in the spatial directions using a Gaussian form, while q(x) 
corresponds to a local field. Additionally, the subscript + on the quark fields implies that they are 
explicitly projected onto positive energy modes, that is u!|_(x, t) = (1 + 74 ) 14 * (x, t). The superscript 
(P, S ) on an interpolating operator (and hence the correlation function) indicates a point or smeared 
form, respectively, C = *7472 is the charge conjugation matrix, and O p = Op ; 74 . Neutron 
correlation functions are constructed from those of the proton by interchanging u o d. The 
correlation functions with the quantum numbers of nuclei, constructed using the methods discussed 
in detail in Refs. [33i, 55), are built recursively using sink-projected nucleon “blocks” that involve 
either smeared or local fields. For the present calculations, zero momentum states are built from 
zero-momentum blocks, although more complicated constructions can also be considered. 

As suggested in Ref. [56], in order to study the proton in a magnetic field, it would be more 
appropriate to use interpolating operators that project onto the lowest-lying Landau level, rather 
than projecting onto three-momentum plane-waves. This would enhance the overlap of the inter¬ 
polating operator onto the lowest, close-to-Landau energy eigenstate and suppress the overlap onto 
higher states. However, it is unclear how to generalize such a framework to nuclei constructed from 
nucleon blocks. While single-hadron blocks provide a good basis for the construction of correlation 
functions of nuclei in the absence of background fields, this will not necessarily be the case in a 
magnetic field. The problem is exemplified by 3 He, a j = \ nucleus comprised of two protons and 
a neutron. Assuming a compact state (which it has been shown to be at the quark masses used in 
this work through calculations in multiple lattice volumes [33]), the wavefunction of the nucleus is a 
Landau level determined by its electric charge of Q = 2e, while that of the proton is a Landau level 
determined by Q = e, and that of the neutron is a momentum eigenstate. Proton blocks could be 
constructed by projection onto a given Landau level at the sink, leading to presumably improved 
signals for the proton correlation functions. However, combining Landau-projected proton blocks 
and momentum-projected neutron blocks will not necessarily produce interpolating operators that 
couple well to the 3 He Landau levels. For larger nuclei, this problem becomes even more severe. 
There are interesting directions that could be pursued in this regard, for example constructing 
single hadron blocks that are projected onto the Landau levels of the “target” nucleus. However, 
a priori it is difficult to ascertain how efficacious such approaches will be. In the current work 
the same interpolating operators used for zero-field calculations are used for calculations in the 
presence of background magnetic fields. Although these are not ideal interpolating operators, they 
are not orthogonal to the expected eigenstates, and the results extracted from this naive approach 
serve as a benchmark for more sophisticated methods that can be explored in future investigations. 

The correlation functions investigated in this work with the quantum numbers of nuclei are of 
the form 


C h , jz (f,B) = <0| f h , jz \ol% p \t-0)} f h , 




|0)i 


( 6 ) 


The exact source and sink interpolating functions, fh-j z and fh-j z , depend on the quantum numbers 
of the nucleus, and are defined implicitly through the recursive procedures of Ref. m • For nuclei 
with non-zero total spin j, the z-component of spin, j z , is made explicit as each magnetic sub-state 
is studied individually. On each QCD gauge configuration, up- and down-quark propagators are 
generated for each of the seven magnetic field strengths from 48 uniformly distributed Gaussian- 
smeared sources. The position of the first source was randomly chosen and the remaining sources 
were placed on a regular, three-dimensional grid relative to that location. The sources locations 
were selected after the background magnetic field was applied and hence there was no correlation 
between the source location and the position of the zero of the vector potential. As the calcula¬ 
tions from the different source locations on each gauge field are averaged over, the dependence of 
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Ch;j z (t', B) on the source location Xj is suppressed. This location averaging effectively projects the 
source interpolating operator onto zero momentum and is discussed in detail in Appendix [Aj In 
most cases, two correlation functions are constructed for each nuclear state using the smeared and 
point sink interpolators, although for larger nuclei there are more possibilities than are calculated. 


C. Magnetic Field Strength Dependence of Energies 


In a uniform background magnetic field, the energy eigenvalues of a hadron, h, either a nucleon 
or nucleus, with spin j < 1 polarized in the z-direction, with magnetic quantum number j z , are of 
the form 


E h . jz (B) = y/M% + Pl + (2n L + l)\Q h eB\ - \x h B - 2^ M 0 ) |B | 2 - 2^^ {2 \f ij B l B 3 ) + ... , 

(7) 


where is the mass of the hadron, Tj is its momentum parallel to the magnetic field, Qh is 
its charge in units of e, and ul is the quantum number of the Landau level that it occupies. 
For a nucleon or nucleus with spin j > there is a contribution from the magnetic moment, 
Hh, that is linear in the magnetic held. The magnetic polarizability is conveniently decomposed 
into multipoles, with j3 h = P^ 1 ^ denoting the scalar polarizability and P^ 1 ^ denoting the tensor 
polarizability (the latter contributes for hadrons with j > 1 ). Tjj is a traceless symmetric tensor 
operator which, when written in terms of angular momentum generators, is of the form 


T ■ — 


a a a a zS a 2 

JiJj + JjJi ~ Oij J 


( 8 ) 


and (...) in Eq. ([t]) denotes its expectation value . 4 Note that the polarizabilities defined here 
represent the full quadratic response to the held and differ from other conventions used in the 
literature where Born terms are explicitly removed (for a discussion, see e.g. Ref. flT?]). The ellipses 
denote contributions that involve three or more powers of the magnetic held. The spin-averaged 
energy eigenvalues project onto the scalar contributions, 


(E h ( B)> = —H £ E h , jz ( B) = ^Ml + Pl + {2n L + l)\Q h eV\ - 2vr/3^ 0 ) |B | 2 + ... ,(9) 

jz = j 

where the ellipsis denotes contributions of C9(|B | 4 ) and higher. For spin-j states, the energy 
difference between j z = ±j isolates the magnetic moment at lowest order in the expansion. Other 
combinations of the energy eigenvalues of the individual spin components can be formed to isolate 
higher multipoles. 


III. RESULTS 

A. Extraction of Energy Levels 

With the background magnetic held given in Eq. ([2]), well-defined energy levels exist for each 
value of the magnetic held strength. In order to determine the magnetic polarizabilities, energy 


4 For a magnetic field aligned in the z-direction, it follows that ( TijBiBj) = {T ZZ B 2 ) = (j? — | j(j + 1)) B 2 . This 
vanishes for both the j = 0 and j = \ states, and takes the values ( TijBiBj) = | B 2 for j = 1, j z = ±1 states and 
(TijBiBj) = — | B 2 for j = 1 ,j z =0 states. 
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eigenvalues are determined from the appropriate correlation functions, the Ch-j z {t: B) defined in 
Eq. ([6]). The individual correlation functions associated with each state in each magnetic field are 
examined, and the time intervals over which they are consistent with single exponential behavior 
are determined. Effective mass plots (EMPs) associated with representative correlation functions 
obtained in the magnetic fields with n = 0,1,—2,3 are shown in Fig. [TJ Having identified these 
time intervals, the main analysis focuses on ratios of these correlation functions constructed as 


Rh;j z {t',B) — 


Ch-,j z {t\ B = 0) 


t—> oo 


Zh-'i 


(B) e ~ SEh »* 


(B )t 


( 10 ) 


where SEf l - l j z { B) = Eh-j z { B) — E/ i; j z (0) is the energy difference induced by the magnetic held, 
and Z} v j z (B) is a time-independent, but field-dependent, overlap factor. To be specific, h = 
p, n, nn, d, pp, 3 He, 3 H and 4 He are considered in all magnetic substates. It is advantageous 
to work with ratios of correlation functions because fluctuations present in the energies extracted 
from individual correlation functions cancel to a significant degree in the ratio. 3 The energy shifts 
are extracted from these ratios in the time regions in which the contributing individual correlation 
functions show single state dominance by either directly fitting the ratio or, alternatively, by fitting 
the effective mass resulting from the double ratio 




1 

tj 


log 


f E h-,j z {t\ B) \ 
\ R h-,j z (t + tj;B)) 


t—¥ OO 


SE h;jz ( B), 


( 11 ) 


where tj represents a temporal offset that can be chosen to optimize energy extractions. The 
fits are performed using the correlated ^-minimization procedure, with the covariance matrix 
determined using jackknife or bootstrap resampling. A systematic fitting uncertainty is estimated 
by performing fits over multiple fitting intervals within the region of single-exponential dominance 
for a given system. The primary analysis in this work is based on fitting the ratios of correlation 
functions obtained from binning source-averaged measurements into IVj, = 100 blocked samples 
and generating Nb = 200 bootstrap samples from these blocked samples. Alternate analyses are 
also undertaken in which differences include varying the statistical procedures and also performing 
constant fits to effective mass functions formed from different values of tj and other possible 
variations. Consistent results are obtained. 

In order to extract the magnetic polarizabilities, ratios of the correlation functions associated 
with the maximally-stretched spin states can be formed such that the leading magnetic-moment 
contributions cancel, 


Rh{t\ B) 


C h - jz=j (t- B) C h . jz= -j(t] B) 

Ch-,j z =j{t'iQ) Ch-j z =—j(t', 0) 




( 12 ) 


where <5(-E/i(B)) = {Eh( B)) — Mh is the spin-averaged energy shift. Similarly, the spin difference 
between maximal j z = states eliminates the spin-independent terms, leaving only the magnetic- 
moment contribution and 0(|B| 3 ) and higher terms. This has been used to extract magnetic 
moments in Refs. IB 03| using the ratio 


Ai4(t;B) 


Ch;j z =j{t; B ) C h -j z =-j(t-, 0) e —(E h .j(B)—Eh i -j{B))t ' 

Ch;j z =j(t'iO) Ch;j z =-j{t' 1 B) 


(13) 


In the present work, the individual spin states are used to extract the magnetic moments and 
polarizabilities in a coupled fit as the latter quantities are the primary target of this study. However, 
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Note that such ratios are formed after averaging over an ensemble of measurements. 
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FIG. 1: EMPs obtained from the smeared-smeared nucleon and nuclear correlation functions for the 
lowest-four magnetic field strengths. The j z = ±j states are shown for each hadron. 
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State 

a SE h . jz (n) 

h 

jz 

n = 1 

fi = —2 

h = 3 

h = 4 

n = — 6 

h = 12 

p 

1 

2 

0.0032(11)(17) 

0.0839(24)(0) 

—0.0324(22)(19) 

—0.0581(26)(12) 

0.1288(51)(65) 

—0.2495(17)(13) 

p 

1 

2 

0.05372(63)(68) 

—0.0073(16)(7) 

0.1035(34)(23) 

0.1087(37)(42) 

—0.1045(57)(41) 

0.142(4)(16) 

n 

1 

2 

0.01297(32)(19) 

—0.03741(12)(8) 

0.0249(12)(19) 

0.0184(21)(12) 

—0.12694(35)(27) 

—0.02318(91)(64) 

n 

1 

2 

—0.01711(7)(17) 

0.01749(75)(46) 

—0.0584(10)(11) 

—0.0831(17)(6) 

—0.0027(23)(24) 

—0.24212(63) (25) 

nn 

0 

—0.00321(15)(32) 

—0.0146(8)(12) 

—0.0285(31)(21) 

—0.0488(32)(36) 

—0.1147(24)(99) 

—0.2793(27)(26) 

d 

1 

0.0190(16)(74) 

0.0588(34) (44) 

—0.0009(54)(31) 

—0.0262(72)(61) 

0.033(7)(15) 

—0.337(22)(19) 

d 

-1 

0.0398(8)(33) 

0.0169(59)(79) 

0.0523(58)(62) 

0.041(8)(13) 

-0.039(47) (28) 

—0.114(7)(20) 

pp 

0 

0.0490(19)(82) 

0.0679(36)(56) 

0.0536(55)(88) 

0.062(5)(12) 

0.001(11)(34) 

—0.114(5)(15) 

3 He 

1 

2 

0.067(3)(24) 

0.0408(38)(53) 

0.123(8)(10) 

0.126(8)(22) 

—0.028(7)(27) 

0.01(1)(69) 

3 He 

1 

2 

0.034(3) (16) 

0.112(4)(13) 

0.023(7)(16) 

0.0045(85)(94) 

0.112(15)(77) 

—0.262(11)(89) 

3 H 

1 

2 

0.007(1)(16) 

0.100(4)(19) 

—0.027(7)(25) 

—0.058(7)(39) 

0.06(2)(11) 

—0.31(8)(12) 

3 H 

1 

2 

0.058(2)(27) 

-0.0025(35) (67) 

0.125(7)(29) 

0.138(8)(38) 

—0.152(11)(35) 

—0.02(2)(70) 

4 He 

0 

0.056(4)(69) 

0.086(7) (30) 

0.090(14)(35) 

0.098(15)(74) 

0.09(1)(11) 

—0.10(3)(62) 


TABLE I: Extracted energy shifts of the nucleons and nuclei, where the first uncertainty is statistical and 
the second corresponds to the systematic uncertainty obtained from variation of fitting ranges. The j z = 0 
deuteron state is studied separately. 


the magnetic moments have also been extracted from the spin-difference ratios, Eq. (13) [I], and 
that approach is found to lead to consistent, but more precise results and remains the preferred 
method for extracting the magnetic moments. Given the more complicated nature of the fits 
we perform here in order to obtain sensitivity to the polarizabilities, it is unsurprising that the 
uncertainties on the lower order terms are larger. In what follows, we present the magnetic moments 
that result from the coupled fits for completeness, but use the previous fits to spin-difference ratios 
as the best estimates of these quantities. 

Figures [2] and [3] show the ratios of correlation functions used in the extraction of energy shifts 
for each magnetic field strength and spin component. Results from all six nonzero magnetic field 
strengths are shown. In each case, the associated single exponential fit to the ratio of correlation 
functions is shown, along with the associated statistical uncertainty. Fits over all time ranges in 
[t min ■ t mav ] are considered, where t max = 24 and t m i n is set by requiring consistency with single ex¬ 
ponential behavior of the individual correlators that form a given ratio. The central fit is identified 
as that over the time range with f max — t m ; n > 3 with the lowest correlated x 2 /d.o.f. The standard 
deviation of all fits over subsets of time ranges in the interval [t m i n — 1 , t max ] that have a x 2 within 
one unit of the minimum y 2 (x 2 -A Xmin + 1) taken as the fitting systematic uncertainty. The 
extracted energy shifts are tabulated in Table [T} 

The correlation functions associated with the nucleons and nuclei are highly correlated, and 
therefore differences between the energies of two given states can be more accurately determined 
than those of each state individually. Of particular interest is the difference between the magnetic 
properties of a nucleus and that of its constituent nucleons. To highlight these differences, further 
ratios are constructed, 




El Rh(t, B) 
h&A 


t—> OO 


5Z_ 4 (B) exp 


f <5FU(B) - ^ 5Eh(B) j t 
\ h&A / . 


(14) 


where the nucleus A contains a set of nucleons, h, and the spin indices have been suppressed for 
brevity. Figure [4] shows this ratio of correlation functions for the nn, j z = ±1 deuteron, and pp 
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FIG. 2: The ratio of correlation functions associated with the p , n, nn and pp systems. Results are shown 
for all six field strengths for the smeared-smeared correlators and for both \j z \ = j states for states with 
j > 0. The shaded bands correspond to the statistical uncertainties of the given fit. 


systems. 

As discussed previously, the momentum-projected interpolating operators are not expected to 
provide particularly good overlap onto the low-energy eigenstates of the proton and charged nuclei 
in magnetic fields, which are expected to more closely resemble Landau wavefunctions. Indeed, 
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FIG. 3: The ratio of correlation functions associated with the d, 3 He, 3 H and 4 He states. Results are shown 
for all six field strengths for the smeared-smeared correlators and for both \j z \ = j states for states with 
j > 0. The shaded bands correspond to the statistical uncertainties of the given fit. 


the interpolating operators are found to overlap most strongly with states other than the lowest 
Landau level, as will be discussed below and in detail in Appendix [Aj In Fig. [5j the ratio of overlap 
factors of the extracted state at nonzero and zero background field strengths are shown. For the 
neutral states, the overlap is only weakly dependent on the field strength, but for charged states, 
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FIG. 4: The ratios of correlation functions defined in Eq. (14) for the nn, d{j z = ±1) and pp systems. 
Results are shown for all six magnetic fields for the smeared-smeared correlation functions. The shaded 
bands correspond to the statistical uncertainties of the given fit. The deuteron spin states are averaged for 
simplicity. 


the overlap rapidly decreases with increasing magnetic field strength. This indicates that care must 
be taken in interpreting the extracted states. It does not mean that unrelated states (those that 
are not continuously related as a function of the magnetic field) are being extracted for different 
field strengths, but instead the overlap onto the given state is decreasing. 6 


6 Similar effects have been seen in constructing multi-pion correlation functions from combinations of pion in¬ 
terpolating operators built with differing momenta. Exponentially smaller overlaps were observed, although a 
consistent energy could be extracted i58l f59i . In background electric fields, the overlap factors of momentum 
projected interpolating operators were also seen to decrease significantly with the applied field strength mum. 
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FIG. 5: The magnetic field strength dependence of the ground state overlap factors, defined in Eq. (10). 
For j z ^ 0 states, both spin states are shown as the different colored points offset slightly for clarity. 


B. The Coupled j z = I z = 0 Two-Nucleon Channel 

The j z = I z = 0 channel is special in that the presence of the magnetic field breaks isospin 
symmetry through the charge matrix and also introduces explicit spin dependence to the low- 
energy effective Hamiltonian. Consequently, the background magnetic field couples the j z = 0 
deuteron and I z = 0 dinucleon states, and a more complicated analysis is required to extract the 
essential physics. The energy eigenvalues of this system result from diagonalizing a 2 x 2 truncated 
Hamiltonian in the basis formed from the ’So and the :i S\ states. 7 For non-vanishing magnetic 
fields, the off-diagonal elements of this truncated Hamiltonian receive contributions from magnetic 
transitions between the 3 Si and ’So channels induced by the nucleon isovector magnetic moment 
and short-distance two-nucleon interactions with the magnetic field resulting from meson-exchange 
currents in the context of potential models or local two-nucleon current couplings in effective held 
theory (EFT). As the nucleon isovector magnetic moment is large, the energy splittings between 
these states are significant for the magnetic held strengths employed. 

The pionless EFT (EFT(^f)) can be used to describe the low-energy strong and electroweak 
interactions of two-nucleon, three-nucleon and multi-nucleon systems [62] , 1 63] . It provides a sys- 


' For the small magnetic field strengths considered in this work, the gap to excitations is significant, and such 
excitations can be neglected. Additionally, due to the tensor interaction, the j — 1 state involves both 5-wave and 
D -wave contributions. The 3 Si- 3 Di coupled channels are truncated down to the 3 5 i channel because the deuteron 
and dineutron are close in energy in the absence of a magnetic field, and the deuteron is predominantly 5-wave 
(at least at the physical pion mass). For a more detailed discussion of the deuteron in a FV, see Refs. [60l [61j . 
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tematic way to include the gauge-invariant electroweak interactions that are not related to strong- 
interaction S-matrix elements through local multi-nucleon operators. While conventionally formu¬ 
lated in terms of four-nucleon interactions with insertions of derivatives (as well as higher body 
interactions), EFT(^) can be fruitfully defined in terms of dibaryon fields, permitting a dramatic 
simplification in calculations beyond leading order (LO) in the expansion [64]. The Lagrange 
density describing the two-nucleon electromagnetic interactions at LO and next-to-leading order 
(NLO) in EFT(/) using dibaryon fields is [ 6 i] 


C = 


2 Mm 


lY f [k 0 + Kir 3 ] S • BN - 




kq ' ^ tj 1-3k 


+ 


^ + h.c. 


M n ^/rir 3 L 0 


(15) 


where ti are the SO(3) vector components of the 3 S'i dibaryon field and s 3 is the I z = 0 , 1 So 
dibaryon field, Mn is the nucleon mass and El is the spin operator. The effective ranges in the 
[So and 3 S\ channels are denoted as ri and r 3 , respectively while Ko and k\ are the isoscalar 
and isovector magnetic moments of the nucleon. The NLO interactions, described by dibaryon 


operators coupled to the magnetic field, are accompanied by the coefficients Z 12 in Eq. (15). 

In Ref. [65], it was recognized that LQCD calculations employing background magnetic fields 
could be used to extract the deuteron magnetic moment, and the rate for low-energy np —> d-y 
radiative capture, by determining the energy eigenvalues of the two-nucleon systems {65 s 66 ]. The 


deuteron magnetic moment is extracted from the energies of the j z = ±1 states in the background 
fields, while the np —> d'y radiative capture cross-section is determined from the nucleon isovector 
magnetic moment and the value of l\ determined from the energies of the two j z = 0 states in the 
coupled — 3 <Si np sector. This latter combination is probed through the determinant condition 


p cot <5i — 


S+ + 5- 

2irL 


p cot 63 — 


S+ + S _ 
2irL 


\eB\h S+ - S- 


1 2 


27 tL 


(16) 


where 5i i3 are the phase-shifts in the x 5o and 3 5i channels, respectively. Solutions to this equation 
correspond to the energy eigenvalues of the system, with the functions S± given by 


S ± = S[^(p 2 ± \eB\ Kl 


(17) 


where 


I n l < A 


S(v) = E 


— 47tA 


n^o 


|n| z — p 


(18) 


is the three-dimensional Riemann-zeta function associated with the irreducible representation 
of the cubic group [671169] . 

At the quark masses used in these calculations, the deuteron and bound dineutron are ap¬ 
proximately degenerate [33], and have scattering lengths, ai )3 , and effective ranges, ri i3 , that are 


8 In this expression, I 2 has been replaced by I 2 — to make explicit the deviation of the deuteron magnetic 
moment, pd = ^ + 1 ^ 7 ° r 3 l) i f rom th e single nucleon contribution. 
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numerically close to each other (ai ~ <23 = a and 77 ~ r 3 = r) [34] and hence ~ 83 = 5 . 9 Because 

|eB|Zi 


of this, Eq. (16) simplifies to 


p cot 5 = —- S± ± 
ttL 


(19) 


where both signs should be taken together for the two solutions. Expanding this for small |eB|, 
the shifts in the energies of the two eigenstates are 


I fp |—^ I _ I p | ^ I 

= -F^d ( K i + joh) + ••• = =F (ki + Li) —+ ... 


( 20 ) 


M ' ' M 

where Z f / = 1/-y/1 — y 0 r is the square-root of the residue of the deuteron propagator at the pole and 


the ellipsis denotes terms that are higher order in the strength of the magnetic field. In Eq. (20), 
the deviations of the energy shifts from their naive single-particle values are defined using 


Li = 70 Zj (h +7-K1) . 


( 21 ) 


To numerically study this coupled system, it proves useful to first construct the correlation 
matrix 


C(t;B) = 


C' 3 5 i, 3 5 i(^ B ) C 3 Sl i So (t-,B) 
Ciso , 3 si(*; B ) c'ts’o.ts’o(*; B ) 


( 22 ) 


where the matrix elements C' J 4 ,s(t;B) are generated from source and sink operators associated 
with the A, B £ { 1 <S'o, 3 <Si} channels (which are orthogonal in the absence of the magnetic field). 
The generalized eigenvalue problem, defined by this correlation matrix, can be solved to extract 
the (diagonalized) principal correlation functions [7U], energies and energy differences. That is, 
solutions of the system 


[C(t 0 ; B )] _ 1 / 2 C(t; B)[C(t 0 ; B)]^ = A(t; B)v 


(23) 


are sought, where the eigenvalues are the principal correlation functions A±(i;B) = exp[— (E ± 
AEa^j^)!] with average energy E and energy difference A E 3 s 1 is 0 - The parameter to can be 
chosen to stabilize the extraction but has little numerical effect in the current results. To extract 
the response to a background magnetic field, the ratio of the principle correlation functions 


-^ISiASo(^i B ) — 


A+(t;B) 

A-(t;B) 


t —^OO 


Zexp [2 AE 3 Slj i So t] , 


(24) 


permits a refined determination of the energy difference A Ez Sl 15 , significantly reducing correlated 
fluctuations, where Z is a t-independent constant. 

Figure 6 ] shows the EMPs of the original correlation functions of the coupled-channel system 
in Eq. (|22) according to their source and sink type. This figure also shows the EMPs constructed 
from the principal correlation functions that are determined by solving the generalized eigenvalue 
problem, Eq. (23), for to = 5. The diagonalization of the matrix of correlation functions in 


Eq. (22) is particularly effective in this case because the states are orthogonal in the limit of 


vanishing magnetic field. In most cases, plateau behavior is visible in both principal correlation 


9 The difference in binding energies is A3£ 1: ig 0 = Ei So — E-j s ^ = 5.8(1.4) MeV uS]; provided the difference in energies 
is small compared to the shifts induced by the magnetic field, it can be neglected. If it cannot be neglected, the 
determinant condition must be solved numerically. 
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FIG. 6: Results from nucleon-nucleon smeared-smeared correlation functions in the mixed j z = I z = 0 
sector. The left panels show the effective masses of the elements of the 2x2 matrix of correlation functions, 
with each quartet of plots corresponding to a different magnetic field strength. In the right panels, the EMPs 
of the principal correlation functions resulting from solving the associated generalized eigenvalue problem 
are shown. 


functions, indicating that the lowest two eigenvalues of the system can be extracted. Given this, 
focus is placed on the ratios R3g 1 ig 0 (t;B) in the region where the principal correlation functions 
are consistent with single exponential behavior. Figure [7] shows the ratios for all magnetic field 
strengths along with the associated single exponential fits. Analysis of these ratios in the coupled 
system is performed with the same methods used to analyze the ratios in the unmixed channels. 

As in Eq. (14), the calculated correlation functions associated with nucleons and nuclei share 
the same quantum fluctuations, to a large degree. This makes it possible to determine differences 
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FIG. 7: The ratio of correlation functions, Rag ± 13 (L B) determined from the principal correlation functions 
for to = 5 for all six magnetic field strengths used in this work. Fits to the correlation functions are also 
shown with uncertainties represented by the shaded region. 


between properties of the np system and those of a free neutron and proton with more precision 
than the individual properties. In the current context, the ratio 


B ) 


B ) 

AR p {t-B)/AR n (t-B) ’ 


(25) 


decays with a characteristic exponent 2 AE3g 1 15 (B) — (E p ^ — E p ±) + ( E n y — E n = 2\eB\L\/M + 
0 (|B| 3 ), permitting direct access to deviations from single nucleon physics, where the ; B) 

are given in Eq. ( |12[ ). Figure [ 8 ] shows these ratios for each field strength, from which the energy 
shifts can be extracted with remarkable precision. 



0 5 10 15 20 0 5 10 15 20 0 5 10 15 20 

t/a t/a t/a 


FIG. 8: The ratio <5i?35 1 is 0 (t; B) computed from the principal correlation functions with to = 5, divided 
by the appropriate isovector combination of the spin differences of the single nucleon correlation functions 
are shown for all six magnetic field strengths used in this work. Fits to these ratios, and the associated 
uncertainties (the bands), are also shown. 


C. Magnetic Field Strength Dependence: General Strategies 

Having extracted the energies and energy differences as a function of the magnetic field strength, 
the remaining task is to use them to determine the magnetic properties of the nucleons and nuclei 
through fits to the expected forms shown in Eq. 0 - The fits and extracted properties of each 
nucleon and nucleus are presented individually in the following subsection. Here, the general 
features of the analysis, and the difficulties encountered in confronting Landau levels, are discussed. 
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In dimensionless units, the form used for the fits to the ground states (B = Be z and T) 

2e 


= 0 ) is 


a SE h . jz = \Ja 2 M 2 + (2 n L + l)Q h a 2 \e Bj - aM h - 


ClMjy 


fi h j z a \e B 


2t r 


a?M 2 N {M A -M N ) ‘ 

+jzlh (a 2 \e Bj ) 3 + 5 h (a 2 |e Bj ) 4 , 


Ph + Ph \j z ~ + 1)) 


(a 2 |e B|) S 


(26) 


where the fit parameters are 


til, 


Ph = 


Hh = Vh- 


Mn_ 
2e : 


7 h, 5 h 


M 2 (M A — Mat) ^(Afo) 

Ph 


2 (2) _ — Mn) a (M2) 

Ph ~ 72 Ph 


(27) 


and a 2 |e B| = is the dimensionless field strength. This extends to higher orders in |e Bj than 
the form given in Eq. ([t]) , providing for estimates of fitting systematic uncertainties in the extraction 
of the magnetic moments and polarizabilities resulting from the choice of fit form. The hadron 
masses, aMh , are taken from our previous studies on this ensemble of gauge-field configurations [33 ) 
and are known precisely (aM& = 1.3321(10)(19) on this ensemble). For uncharged states, ul does 
not enter the fit, and for states with j = 0, the parameters p, p y h ; and 7 are absent. As in 
Ref. [lj, the extracted magnetic moments are expressed in terms of “natural nuclear magnetons” 
( 11 NM) defined with respect to the nucleon mass at the quark masses used in the calculation. The 
polarizabilities are given in terms of the natural dimensionless polarizability e 2 — M/v) 
(given the expected dominance of the A-resonance, this is the appropriate scale for the magnetic 
polarizabilities), but are also presented in physical units in the conclusion. A physical interpretation 
of the higher order parameters is not provided, and they are used only to control the systematic 
uncertainties in the magnetic moments and polarizabilities extracted from the fit. 

In performing fits, the same bootstrap sets of extractions of the energy shifts are used at each 
magnetic field strength in order to exploit the correlations between them. The ensemble averages 
of the energy shifts are used to obtain the central values of fit parameters describing the magnetic 
field strength dependence. An ensemble of fits to the bootstrap data set is used to obtain the 
associated statistical uncertainties. To propagate the systematic uncertainties from the fits to the 
ratios of correlation functions into the field dependence analysis, the bootstrap sets of energy shifts 
are spread away from their mean by the ratio of the quadrature-combined statistical and systematic 

uncertainties to the statistical uncertainty. That is, for an energy variable E with an ensemble of 

_ n b 

extracted bootstrap values {£)}, mean value E = jr- Eb, and statistical and fitting systematic 

B b =1 

uncertainties 6 E Sta t and SE sys , the spreaded bootstrap ensemble values {£)} are taken to be 


j5E 2 tat + 5E 2 ys 

Ei = E+ A-_- (Ei - E ), 


5E. 


stat 


(28) 


and it is these quantities that are used in subsequent analyses. The highly correlated nature of 
the results obtained at different field strengths (that is, results obtained at different magnetic field 
strengths on a given configuration have correlated statistical fluctuations) makes this approach 
essential. 

An important aspect of the present analysis is to address the range of magnetic field strengths 
for which the fit forms in Eq. (26) describe the energy shifts. This is addressed by varying the 
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number of magnetic field strengths used in the fits, fitting results obtained at |h| < h max with 
h max = 2, 3,4, 6,12. The complexity of the fit form is also varied by either using the full form or 
by setting 5 = 0, 5 = 7 = 0 or <5 = 7 = /? = 0. Further, either the smeared-sink or point-sink 
correlation functions are selected for fitting. For j = 0 states, a total of 4x2x2 = 16 different 
fits are considered for each bootstrap ensemble. For the j > 0 states, coupled fits were performed 
to the magnetic field strength dependence of the j z = ±j spin states and 5 x 4 x 2 = 40 different 
fits are considered for each bootstrap ensemble. A large number of fits successfully described the 
results with an acceptable y 2 /d.o.f., although some did not. The central values and uncertainties 
in the extracted parameters are evaluated from the distribution of the results for the acceptable 
fits, and are taken as the 50 th quantile and the 17 th -83 rd range of quantiles, respectively. As an 
additional check, we have used the Bayes information criterion for a given fit to assess those that 
are acceptable and find uncertainties that are consistent with those defined from the y 2 . 

A further complication arises from the Landau-level nature of the eigenstates and the sub- 
optimal projection of the interpolating operators onto them. In every case, correlation functions 
that have single exponential behavior over a significant range of time-slices are found. However, 
it is clear that these states do not correspond to the lowest Landau level. Expanding the field 
dependence in the non-relativistic limit, the magnetic moment contribution cancels in the spin- 
averaged energy shift, but a linear contribution survives with a coefficient determined by ul, 

a [E h , jz { B) + E h ._ jz ( B)] = + 1) + 0(|B| 2 ), (29) 


where the masses are precisely known from previous studies |33j . From examining the small field 
shifts for charged states, it is found that til ^ 0 in all cases. Thus, the interpolating operators 
overlap strongly onto excited states of the system and presumably will relax to the ground state only 
at large Euclidean times. Because of this, til is treated as a fit parameter and the fits themselves 
are used to determine which Landau levels are dominating the various correlation functions. In 
the limit of vanishing lattice spacing, and neglecting structure effects, the allowed values of til are 
integers and a somewhat complicated approach to fitting is required. Two alternate procedures 
are considered. In the first approach, the lowest magnetic field strengths are used to determine 
the linear term in the field-strength dependence, which is used to identify the integer value of 
til that is most consistent with the numerical results. This value is held fixed and then used in 
further fits utilizing the full form of Eq. (26). In the second approach, rip is first treated as a 
real-valued fit parameter and the full fits are performed. Then, after considering the different fits, 
the integer n l closest to the mean of successful fits is chosen and held fixed in the final set of fits . 10 
An additional systematic is assessed by combining sets of fits (varying fit forms, data ranges, and 
types of correlation functions) with hp —> hp ± 1 into the full suite of fits (for charged, j > 0 
states, a total of 120 different fits are considered for each bootstrap ensemble). Both choices of 
Landau-level procedures lead to consistent results after these systematic uncertainties are taken 
into account. 

A related systematic uncertainty that is considered is the potential ambiguity in identifying the 
Landau level associated with the plateau at each magnetic field strength. In the fitting forms, it is 
implicit that the energies of a nucleon or nucleus result from a single Landau level for the range of 
magnetic fields that are considered. This is expected to be the case at small magnetic fields, but 
might not be valid at larger magnetic fields. To explore this issue, fit forms with different hp for 
different magnetic fields are considered. Keeping only the results from the lowest four magnetic 


10 


The extractions of magnetic moments and the np —r dy cross-section are independent of these complications as 
the Landau-level contributions cancel in the energy differences between spin states. 
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fields, and allowing different values of til for each magnetic field does not result in acceptable fits 
except when the h^’s are all the same. This leads to confidence in the assumption that the same 
Landau state is providing the energies that are dominating the fits. 


D. The Magnetic Properties of Nucleons and Nuclei 

1. The Neutron 

The neutron correlation functions and their ratios for each spin component and magnetic held 
strength used in this analysis, along with the associated fits, are shown in Figs. [l]and[2j and the 
energy shifts extracted from these functions are presented in Table [T| Figure [9] shows the energy 
shifts of a spin-up and spin-down neutron as a function of the background magnetic held strength. 
The two spin states behave quite differently in the presence of the magnetic held. The energy shift 
of the spin-down state (negatively shifted as the neutron magnetic moment is negative) responds 
almost linearly to the magnetic held, even out to |eB| ~ 0.71 GeV 2 (n = 12), while the response of 
the spin-up state exhibits signihcant nonlinearities even for modest magnetic helds. The behavior 
of the spin-up state is reminiscent of the lower level in a two-state system with an avoided level 
crossing. Given the expected tower of QCD excitations of the nucleon, the observed behavior of the 
spin-up state is consistent with the magnetic held inducing mixing between the spin-up neutron and 
higher-lying states. Such mixing is expected from quark-hadron duality, and LQCD calculations 
that also probe the response of excited states to the magnetic held could be used to investigate 
this further. This behavior implies that spin-dependent polarizabilities are highly correlated with 
spin-independent polarizabilities, and it will be interesting to learn if this pattern persists as the 
quark masses are brought closer to their physical values. 

As discussed in the previous subsection, a large number of different hts involving varying ranges 
of held strengths and with a variety of functional forms are performed in order to the analyze 
the energy shifts and determine the magnetic moment and polarizability. The 17 th and 83 rd 
quantile range of all successful hts to each data range are shown as the shaded regions in Fig. [9} 
Separate (overlapping) regions are shown for each data range n < n max for which successful hts 
were found. Note that the hgure shows the results from smeared-smeared correlation functions only, 
but the hts that are considered also include those involving the energies extracted from smeared- 
point correlators. The hgure also shows the probability density functions (PDFs) generated from 
combining the central values of all successful hts (considering hts over the energy shifts extracted 
from each of the bootstrap ensembles) for the two relevant parameters, jl and (3. For hts involving 
additional parameters, these are integrated over, while for linear hts just involving the magnetic 
moment, these are ignored in determination of the PDF for (3. Analysis of the suite of hts that are 
detailed in the previous subsection yields a neutron magnetic moment and polarizability of 

An = -1.972 (”) (0.059), (30) 

Pn = 0-198 (“) (0.010), (31) 

where the hrst uncertainties combine the statistical and systematic uncertainties from the extrac¬ 
tion of the energy shifts, as well as the systematic uncertainty from the ht to the magnetic held 
strength. The second uncertainty estimates the effects of discretization and finite volume effects; 
as discussed in the conclusion this is assessed to be a 3% multiplicative uncertainty on magnetic 
moments and a 5% multiplicative uncertainty on polarizabilities. The above results are presented 
in the natural dimensionless units, and the values of the magnetic moment and polarizability in 
physical units are subject to additional uncertainties from the lattice scale-setting procedure, which 
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FIG. 9: Results for the energy shifts of a spin-up (upper points) and a spin-down (lower points) neutron 
in a background magnetic field. The central 68 th quantile of successful fits is shown as the shaded bands. 
Different overlapping bands are shown for fits over the different ranges of h. The lower panel shows the 
probability-density functions for the relevant fit parameters /x and /3, with the vertical lines indicating the 
central value and uncertainties. 


are discussed in the conclusion. The magnetic polarizability and magnetic moment of the neutron 
have been calculated previously with LQCD over a range of light-quark masses [38, 23j albeit with 
large uncertainties. The calculated magnetic moment is consistent with previous calculations at 
similar quark masses, and the value of /3 n is also consistent with previous calculations [35], 11 


2. The dineutron 

At these unphysical quark masses, the dineutron (in the ’So channel) is a bound state, with a 
binding energy of B nn = 16(5) MeV [33]. As it is electrically neutral, comprised of two neutrons 
in the 1 S'o channel with positive parity, the dineutron provides the simplest nuclear system with 


n 


The authors of Ref. [43] report difficulties in identifying ground states. 








23 



0.20 0.25 0.30 0.35 0.40 


FIG. 10: The energy shifts of the dineutron as a function of the background magnetic field strength, h. 
The details of the figure are as in Fig. [9] The lower panel shows the PDF for the dineutron polarizability. 


which to explore the effects of binding on magnetic properties. This system is discussed before 
proceeding to states that are electrically charged and therefore complicated by the presence of 
Landau levels. 

Figures [l] and [2] show the dineutron correlation functions and the ratios of correlation functions 
for each field strength, along with fits to the time dependence of the ratios. The energy shifts 
extracted from the ratios of correlation functions are given in Table]]] for each field strength, and 
Figure [T0| shows these shifts. Combining all of the attempted fits to the energy shifts, as described 
in detail for the neutron, yields a magnetic polarizability of 


P nn = 0.296 ( +° ° o ; 9 8 ) (0.015), 


(32) 


where the uncertainties are as for the case of the neutron, and the result is presented in the 


dimensionless natural units of the system, defined in Eqs. (26) and (27). 


This polarizability is significantly smaller than twice the single neutron polarizability with 
5/3 n n = Pnn—ZPn ~ 0.1. This difference can also be obtained from the ratio 5R nn fl(t ; B) in Eq. ( fl4| 
that probes the difference directly in a correlated manner. In the large time limit, the exponential 
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FIG. 11: Results for the difference between the energy shifts of the dineutron and a spin-up and spin-down 
neutron as a function of the background magnetic field strength. The green shaded regions are the result 
of the suite of fits to the field strength dependence. The horizontal red-shaded region shows the breakup 
threshold for the dineutron, above which the ground state of the system would be two neutrons in the 
continuum in the ISo channel. 


decay of this ratio is governed by the energy difference 5E nn ( B) — 6 E n i (B) — 5E n _i (B). These 


ratios are dis play ed in Fig. [4] for the different field strengths and the extracted energy shifts are 


shown in Fig. 11 as a function of the field strength. In turn, 8f3 nn is the coefficient of the quadratic 


term in the field strength dependence of this energy difference. Analyzing these energy shifts using 
the same methods as above, leads to 


5p nn = p nn - 2f3 n = -0.070 ( ) (0.004) . (33) 

As discussed above, the neutron spin-down state is magnetically rigid and remains undeformed 
even at large magnetic fields, while the spin-up state is strongly deformed. For the dineutron, the 
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overall energy is lowered in a magnetic field, driven largely by the spin-down neutron. As this also 
lowers the energy of the spin-up neutron, it has a reduced mixing with other states and, therefore, 
becomes more rigid. From Fig. 11 it is apparent that the binding energy of the dineutron (the 


energy required to separate the spin-up and spin-down neutron) is reduced for h< 8, but is larger 
for n = 12, than at B = 0. If this behavior persists at the physical quark masses, it would indicate 
that it is energetically disfavored for neutron matter (or near-neutron matter) in dense stellar 
objects to spontaneously generate a magnetic field through the formation of dineutron pairs. It 
is also interesting to note that at intermediate field strengths, the dineutron system is nearing a 
Feshbach resonance in which the binding energy is approaching zero and the scattering length is 
approaching infinity. 


3. The proton 


The analysis of the proton in a magnetic field is more complicated than that of the neutron 
and dineutron. As discussed previously, the interpolating operators used in this work project onto 
plane waves in all three spatial directions rather than Landau levels which we expect to be closer 
to the eigenstates of the system, so the quality of the correlation functions for charged systems is 
expected to be significantly worse than that for electrically neutral systems. This is indeed the 
case, as can be seen from the EMPs shown in the first row of Fig. [l] in comparison to the neutron, 
the proton correlation functions are of lower quality with plateaus setting in at later times and with 
significantly larger uncertainties. Further, the presence of Landau levels significantly complicates 
the spectrum of charged system and it is clear that the plateaus that are evident do not correspond 
to the lowest Landau level, as discussed above. The Landau level associated with the plateau is 
identified through systematic analysis of the field-strength dependence, as discussed in Section 

mm 

In Fig. [2j the ratio of correlation functions associated with each spin component and field 
strength, and associated fits, are shown. The energy shifts resulting from these fits are given 
in Table |T] for each magnetic field strength and are shown in Fig. 12 
performed lead to a proton magnetic moment and polarizability of 


The suite of fits that are 


a? = 3 - i7 ( )(°-° 9 )> 


(3 p = 0.83 


+ 0.10 


(0.04) , 


(34) 

(35) 


where the uncertainties are as discussed for the case of the neutron, and the results are presented 
in the dimensionless natural units defined in Eqs. (26) and (27). The Landau level makes a 


contribution to the O (L> 2 ) term that is suppressed by the mass of the proton, and its main 
contribution is to the term linear in B where it can be well constrained by the coupled analysis 
of the two spin states. Consequently, the uncertainty in identifying the correct Landau level of 
the system does not lead to a particularly large uncertainty in the extracted value of its magnetic 
polarizability (although the neutron polarizability is considerably more precise). 

The magnetic polarizability of the proton is found to be considerably larger than that of the 
neutron, (3 p — f3 n = 0.63(10) (4), indicating a significant isovector component at this unphysical pion 
mass. Currently, there are no other LQCD calculations of the proton magnetic polarizability with 
which to compare, however, it can be compared with the experimental value. As quoted previ¬ 
ously, /3p hyh = 3.15(0.35)(0.2)(0.3) x 10 -4 fm 3 |3K8] which corresponds to /3p hyb = 0.116(13)(7)(11) 
in dimensionless units. The physical value results from cancellations between pion-loop (chiral 
physics) and A-pole contributions that are both 0(10 x 10 -4 fm 3 ). Since the pion-loop contri¬ 
bution is strongly suppressed at heavy quark masses and the A-pole contribution depends less 
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FIG. 12: Results for the energy shifts of a spin-up (lower points) and a spin-down (upper points) proton 
in a uniform background magnetic field. The details of the figure are as in Fig. [9] The lower panel shows 
the PDFs for the fit parameters p, and $. 


strongly on mass, the size of the magnetic polarizability determined at the SU(3) point is in line 
with expectations. 


4- The diproton 


The diproton is in the same ’So isotriplet as the dineutron and, neglecting the electroweak 
interactions and the difference in mass between the up- and down-quarks, it would have the same 
properties as the dineutron at zero magnetic held. However, the presence of the background 
magnetic held breaks isospin symmetry through the light-quark electric charges, so the diproton 
magnetic properties are expected to be quite different from the dineutron, even neglecting the issue 
of Landau levels. 

Extracting energy differences from hts to the ratios of correlation functions shown in Fig. [2] 


leads to the results shown in Fig. 13 Fitting the energy shifts, as discussed previously, allows for 
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FIG. 13: Results for the energy shifts of the diproton as a function of the background magnetic field 


strength. The details are as those in Fig. 10 The lower panel shows the PDFs for the fit parameter /3. 


an extraction of the diproton polarizability of 


ftp = 0.84 ( ) (0.04) . 


(36) 


As in the case of the dineutron, the correlated ratios of the diproton and the spin-up and spin-down 
proton correlation functions directly determines the difference of energy splittings. Figure [4] shows 
these ratios, leading to the energy shifts shown in Fig. 14 The figure also shows the envelopes of 


the ensemble of acceptable fits that were performed using polynomials of up to quartic order. 

It is clear from Fig. [14] that the magnetic field strengthens the binding of the diproton by a 
significant amount that rapidly increases until h ~ 3 and then remains constant for larger field 
strengths. This behavior is interesting in the context of the suggestion that at the physical quark 
masses, the diproton can overcome the Coulomb repulsion and form a bound state m in a strong 
enough magnetic field, although this argument requires the system to be near unitarity. However, 
the form of the difference is more complicated in this case than for the dineutron because the 
contributions of Landau levels in the diproton and spin averaged protons may be different. The 
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FIG. 14: Results for the difference between the energy shifts of the diproton and a spin-up and spin-down 
proton as a function of the background magnetic field strength. The red-shaded region corresponds to the 
breakup threshold, above which the ground state of the system would be two protons in the continuum in 
the % channel. 


difference in magnetic polarizabilities is therefore estimated in the naive way, giving 

sew = -°- 82 (! o : 3 2 7 )(°- 04 )’ ( 37 ) 

where the uncertainties of the diproton and proton polarizabilities are combined in quadrature. 


5. The deuteron : j z = ±1 

The deuteron is a bound state in the positive parity :i S\- :i D\ coupled channels. In a background 
magnetic field, while the j z = ±1 states remain isolated in the 3 5 i- 3 Di coupled channels (in infinite 
volume), the j z = 0 state mixes with the positive parity 3 5o isotriplet np channel. Here, the focus 
is on the j z = ±1 states which are used to and extract the magnetic moment and a combination of 
the scalar and tensor polarizabilities. The j z = I z = 0 coupled states are addressed in the following 
subsection. 

Figure [l] shows the effective masses resulting from the j z = ±1 deuteron correlation functions 
and Fig. [3] shows the ratios of these correlation functions, along with fits to their time dependence. 
The energy shifts extracted from these ratios are shown in Fig. |To} Analysis of the field strength 

dependence through a suite of coupled fits to the two spin states, as discussed above, leads to a 

magnetic moment and polarizability of 

fi d = 1.41 ( +0 o 2 2 8 5 ) (0.04), (38) 

Pd M0) + \Pd M2) = 0-70 ( + 0.23 ) (°- 04 ) • (39) 

As the deuteron has j = 1, both the scalar and tensor polarizabilities contribute to the quadratic 
dependence on the magnetic field strength, as presented in Eq. ([ 7 ]). 12 


12 With further analysis, the G(|eB| 2 ) shifts in the j z = 0 np coupled system should determine an orthogonal 
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FIG. 15: Results for the energy shifts of the deuteron in the j z = ±1 states as a function of the background 
magnetic field strength (the lower points correspond to the j z = +1 state). The details of the figure are as 
in Fig. 10 The lower panel shows the PDFs for the fit parameters fi and (3. 


The sum of the proton and neutron magnetic polarizabilities at this pion mass is f3 p + f3 n ~ 
1.02 ^ ^ (0.05), so the deuteron in the j z = ±1 states is somewhat more magnetically rigid 

than the sum of its constituents. While they cannot be separated from this result alone, the 
nuclear forces and gauge-invariant electromagnetic two-nucleon operators are responsible for this 
difference. If the difference persists at the physical quark masses, this would suggest that the 
extraction of the neutron polarizability from experiments on the deuteron is problematic. Figure 


as a function of the field strength. As in the case of the dineutron, the magnetic field pushes 
the j z = ±1 spin states of the deuteron towards threshold and at h ~ 5, the deuteron becomes 
potentially unbound before rebinding at larger field strengths. The figure also shows the envelopes 


combination of the scalar and tensor polarizabilities, /3^ M0 ^ — as given in Eq. |t]) , but this extraction is not 

pursued in the present study. 


16 shows the splitting between the j z = ±1 spin states of the deuteron and the breakup threshold 
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of the ensemble of acceptable fits that we perform using polynomials of up to quartic order. As for 
the case of the diproton, the presence of Landau levels that may differ between the deuteron and 
proton complicates the analysis of the field strength dependence and we do not report a value of 
Wd- 



FIG. 16: Results for the difference between the energy shifts of the j z = +1 spin states of the deuteron 
and that of a spin-up neutron and spin-up proton as a function of the background magnetic field strength. 
The red-shaded horizontal band corresponds to the breakup threshold, above which the ground state of the 
system would be a proton and a neutron in the ISq continuum. The green shaded regions correspond to the 
envelopes of the fits discussed in the main text. 


6. 3 He 


At the physical quark masses, 3 He can be thought of, to a large degree, as two protons spin- 
paired in the ] Sq channel and a single unpaired 5-wave neutron. The ground state is positive parity 
with spin-half and is an isospin partner with the ground state of the triton, 3 H. A naive shell- 
model prediction is that the magnetic moment of the ground state of 3 He is that of the neutron 
(with the spin-paired protons not contributing) and that the magnetic moment of the ground 
state of the triton is that of the proton (with the spin-paired neutrons not contributing). The 
experimental values of both magnetic moments deviate only slightly from these naive predictions. 
Recent calculations have shown that this feature persists even at heavier quark masses m , in 
particular, at the pion mass employed in the present analysis. 

The EMPs obtained from the 3 He correlation functions in the background magnetic fields are 
shown in Fig. [l] and the ratios of correlation functions for each spin state are shown in Fig. [3j 
along with fits to their time dependence. The quality of these ratios is inferior to those obtained 
in the one-nucleon and two-nucleon sectors, but strong signals are still evident. The energies that 
are extracted from these ratios, are shown in Fig. 17 Analysis of the field strength dependence of 
the two spin states allows the magnetic responses to be determined, leading to a magnetic moment 
and polarizability of 3 He of 


A 3 He 

= -2.28 

{ +0.59 
V -1.04 

I? 

O 

(40) 

+ He 

= 0.85 ( 

+0.34 \ 
-0.32 ) 

(0.04) , 

(41) 
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FIG. 17: Results for the energy shifts of 3 He as a function of the background magnetic field strength along 
with the fit envelopes. The details of the figure are as in Fig. [lOj The lower energy points correspond to the 
j z = — 1 state, while the upper points correspond to j z = +|7The lower panel shows the PDFs for the fit 
parameters /t and /3. 


in natural dimensionless units. Given that the magnetic moment favors the naive shell model 
expectation, /i3 He = Tm one might naively expect the polarizability of 3 He to arise from the 
polarizability of the diproton combined with that of the neutron. Within the uncertainties, that 
expectation is found to hold. 

The uncertainties in the magnetic polarizabilities of 3 He are sufficiently large that statistically 
significant deviations from the contributions from the one-body contribution are not obtained, and 
hence we have no meaningful constraint on the MEC contributions. 


7. The triton 


As in the case of 3 He, the ratios of the triton correlation functions are significantly less well- 
defined than those in the one-body and two-body sectors. The energy shifts extracted from the 
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correlation functions are shown in Fig. 18 


Fits to the magnetic field strength dependence of the 



FIG. 18: Results for the energy shifts of 3 H as a function of the background magnetic field strength, 
along with envelopes of fits. The details of the figure are the same as in Fig. [To] The lower energy points 
correspond to the j z = +\ state, while the upper points correspond to j z = — |~The lower panel shows the 
PDFs for the fit parameters jl and (3. 


energies of the two spin states enable an extraction of the magnetic moment and polarizability of 
the triton of 

Ash = 3.32 (+° 0 ™) (0.10), (42) 

^ = 0.40 ( +“; 2 / 7 ) (0.02) . (43) 

The value of the triton polarizability is considerably smaller than the naive expectation of the sum 
of the polarizability of the dineutron and of the proton, f3 p + /3 nn = 1.12 ^ ^ and this differ¬ 

ence could potentially be used to provide a constraint on two- and three-nucleon electromagnetic 
interactions. 
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FIG. 19: Results for the energy shifts of 4 He as a function of the background magnetic field strength. The 


details of the figure are as in Fig. 10 The lower panel shows the PDF for the fit parameters /?. 


8. 4 He 

The 4 He nucleus has the quantum numbers of two protons and two neutrons in a spin-zero, 
even-parity configuration. The energy of the ground state has been determined at unphysical 
quark masses in previous LQCD calculations [MUZUHJ, and at this pion mass it is bound by 
i?4 He = 107(24) MeV [33j. While it has no magnetic moment, it can be polarized by electromagnetic 
fields. 

The EMPs obtained from 4 He correlation functions in the background magnetic fields are shown 
in Fig. [lj and the ratios of correlation functions are shown in Fig. [3j along with fits to their time 
dependence. The energy shifts extracted from fits to these ratios are given in Table [T] and are 
shown in Fig. 19 Analysis of the magnetic field strength dependence of the 4 He energies enables 


an extraction of the magnetic polarizability, giving 


P*He = 0.54 


+0.32 

-0.31 


(0.03) . 


(44) 
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E. The j z = I z = 0 np states and the np —> d'y transition 

The two energy eigenvalues of the coupled j z = I z = 0 np channels are shown in Fig. [ 6 ] for 
each magnetic field strength. In order to extract the : \S\ - 'So mixing, and hence the short-distance 
two-nucleon (MEC) contribution to np —> d'y, the ratios of correlation functions R,3s 1 is 0 (t\ B) in 
Eq. (24) are constructed, and shown in Fig. [7J along with fits to the time dependence. The energy 
shifts extracted from these ratios are shown in Fig. 20, along with the envelope of the ensemble of 
successful fits, from which the linear coefficient is found to be 


ki + Li = 2.74 


+0.07 

-0.05 


(0.07) nNM 


(45) 


where the result is presented in dimensionless units determined by the natural nuclear magneton 
at this pion mass. Fits of up to quadratic order are considered in this analysis. 

To further isolate the short-distance two-nucleon contribution, the ratios SR:ig ] [s 0 (t] B), defined 
in Eq. (25), are formed. By design, the energy shifts extracted from these ratios (see Fig. [ 8 ]) have 
the form 2|eB|Li/M + 0(|eB| 2 ). These shifts are shown in Fig. 21, and performing polynomial 
fits to the dependence on the magnetic field strength leads to 


Li = 0.207 


+0.020 

- 0.020 


(0.006) nNM 


(46) 


Given the isovector magnetic moment and the short-distance two-nucleon contribution, the 
cross-section for the process np —> d'y can be determined near threshold at ~ 806 MeV. Even 
though the LQCD calculations determines these parameters from mixing between bound states, 
the EFT(^f) framework is valid for low-energy scattering states and can be immediately applied. 
It is conventional to use a multipole expansion to define the cross-section for the radiative-capture 
process np -» d'y at low energies [631 ESI [76], 


a(np —y d'y ) = 


253 


e 2 (7o 2 + | P | 2 ) 

m at7oIpI 


\Xmi | 2 + |-^bi | 2 + 


(47) 


where Xmi is the Ml amplitude and Xei is the El amplitude for the process, 70 is the binding 
momentum of the deuteron and p is the momentum of each incoming nucleon in the center-of- 
mass frame. The ellipsis denotes higher-order amplitudes, suppressed by powers of the photon 
momentum. Following Refs. [64L165] . it is straightforward to show that the amplitudes at NLO, 


with the dibaryon parameterization of Eq. (15), are 


Xei = 
Xmi = 


1 


|p| M A r 7 o 


V 1 - 7CT3 (IpI 2 + 7o) 2 
Z d 


-ir + i r i|p| 2 -*|p| 


K i7o 


L 7o + 


1 1 

70-l- -n p 

a\ 2 


, To 7 ' 

+ T 1 


(48) 


where the quantities appearing in this expression are defined in Section IIIB Near threshold, the 
El amplitude is sub-leading and will be ignored. Inserting the extracted values for k±, Li, the 
binding energy from Ref. [33}, and the scattering lengths and effective ranges from Ref. [31], leads 
to a radiative capture cross-section at the SU(3) symmetric point of 


o(np —y dq)! 


m7r~806 MeV 


i7 ( + -“«) > 


(49) 


for an incident neutron speed of v = 2200 m/s, accurate up to NLO in EFT(^). Because of the 
non-linear nature of the dependence of the cross-section on the inputs, the distribution is extremely 
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FIG. 20: Results for the differences in energy shifts between the two j z = I z = 0 np energy eigenstates as 
a function of the background magnetic field strength. The details of the figure are the same as in Fig. [Toj 
The lower panel shows the PDF for the coefficient of the linear field dependence, Hi + L\. 


non-Gaussian; the central value is reported as the 50 th quantile and the uncertainty bounds as the 
17 th and 83 rd quantiles of the full distribution. Improving on this uncertainty requires significantly 
better determinations of the scattering parameters and the binding momentum. At the physical 
point, the cross-section is known to be a(np —> d'y) = 334.2(0.5) mb [77] at this relative velocity, 
which is significantly larger. The short-distance two-body contribution in the calculated cross- 
section (Eq. (49)) is about 10%, just as in the phenomenological determinations. Accounting 
for the significantly different phase space available at the SU(3) point, and the greatly different 
scattering parameters in both channels. The discrepancy in the cross-section is unsurprising. In 
Ref. [29], this result is combined with an analogous result at m ~ 450 MeV to extrapolate to the 
physical point and postdict a cross-section of cr lqcd (np —> d'y) = 332.4 7 ) mb. 
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FIG. 21: The correlated difference between the differences in energy shifts between the two np j z = 0 energy 
eigenstates and those of the isovector nucleon (see Eq. ( |25|) a s a function of the background magnetic field 
strength. The details of the figure are the same as in Fig. |10| 


IV. SUMMARY 

The magnetic moments and magnetic polarizabilities of the lightest few nuclei have been calcu¬ 
lated at a pion mass of m n ~ 806 MeV using LQCD in the presence of background magnetic fields. 
In addition, by considering the mixing of two-nucleon states with j z = I z = 0, the L\ counterterm 
of EFT(^) that governs short-distance two-nucleon contributions to the radiative-capture process 
np —> d'y has been determined. This has then been used to predict the near threshold capture 
cross-section at this pion mass. The success of the calculations presented in this work, and in 
Refs. demonstrate the feasibility of studying the structure of nuclei directly from QCD and 

open the way to a variety of additional QCD calculations of the structure and interactions of light 
nuclei. 

The LQCD calculations have been performed at a single lattice spacing and in one lattice volume, 
and the lack of continuum and infinite-volume extrapolations introduces systematic uncertainties 
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into the results. The effects of the FV used in this work on the binding energies of the light nuclei 
have been explicitly quantified in previous works [33] and found to be small. It is expected that such 
effects in the moments and polarizabilities are of comparable size. An additional uncertainty of 
e~ lL is assigned to the extracted values of nuclear moments and polarizabilities and for simplicity, 
we conservatively use 7 = 7 d ~ 190 MeV for the binding momentum, leading to a ~ 3% uncertainty. 
For the single nucleons, the expected volume effects of 0(e~ m7rL ) are negligibly small. Calculations 
with multiple lattice spacings have not been performed and this systematic uncertainty remains 
to be quantified. The electromagnetic contributions to the action are perturbatively improved 
as they are included as a background field in the link variables. Therefore, the lattice-spacing 
artifacts are expected to be small, appearing at ^(Aq CD a 2 , o^Aqcdo) ~ 3% for Aqcd = 300 
MeV. To account for these effects in dimensionless quantities, an overall multiplicative systematic 
uncertainty of 3% is assigned to the extracted magnetic moments and L 1 , and an uncertainty of 
5% is assessed on all of the polarizabilities, where more complicated effects that compound these 
uncertainties may be possible. For nuclei, these contributions are small compared to the other 
systematic uncertainties. The main results are presented in terms of dimensionless quantities, 
but in Table m we also convert the polarizabilities to physical units using the lattice spacing 
a = 0 . 110 ( 1 ); since the units of polarizabilities are fm 3 , the scale-setting uncertainty corresponds 
to an additional 3% uncertainty that is added in quadrature. Unfortunately, the calculations of the 
individual polarizabilities are incomplete because of the omission of the disconnected contributions 
(the coupling of the external field to the sea quarks), however empirical evidence (471 1521 [53] 
suggests that the omitted contributions will lead to only small modifications that lie within the 
current uncertainties. Confirming this expectation is left to future work. We stress that the 
magnetic moments and the M 1 transition amplitude for np —> d'y are not afflicted by the absence 
of disconnected diagrams (and nor are isovector differences such as f3 p — (3 n at the SU(3) point). 

The magnetic moments and polarizabilities that have been determined in this work and in 
Ref. [T] are summarized in Tables [IT] and III and are shown in Figs. 22 and 23 (the magnetic 
moments calculated from spin splittings in Ref. [1] are the most precise determinations). The 
electrically neutral systems are found to be by far the most precise because the electrically charged 
systems are defined by Landau levels, which have less than ideal overlap with the interpolating 
operators used to form the correlation functions. 


State 

/i[n NM] 

n 

—1.981(05)(18) 

p 

+3.119(33)(64) 

1—1 

A 

II 

^3 

+1.218(38)(87) 

3 He 

—2.29(03)(12) 

3 H 

+3.56(05) (18) 


TABLE II: The results of our previous calculations of the nucleon and nuclear magnetic moments [l] from 
spin splittings at a pion mass of tn, ~ 806 MeV. The first uncertainty is statistical while the second is the 
complete systematic. As discussed in the text, these values are more precise than those determined from 
the more complex analysis required to extract the polarizabilities. 

These results, while not obtained at the physical values of the light-quark masses, are interest¬ 
ing in their own right and suggest important features of these systems. First, our calculations are 
sufficiently precise to determine that the strong interactions between neutrons are such that when 
placed into a magnetic field, the two-neutron system is more magnetically rigid than the sum of the 
individual neutrons. This is consistent with expectations at the physical quark masses based upon 
phenomenological nuclear interactions m, and indicates that it is energetically disfavored for a 
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TABLE III: The magnetic polarizabilities calculated in this work at a pion mass of m, r ~ 806 MeV. An 
additional 5% uncertainty is associated with each polarizability as an estimate of discretization and finite 
volume effects. For the polarizabilities presented in physical units, an additional scale setting systematic 
uncertainty (3%) is included in quadrature in the second uncertainty. 
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FIG. 22: A summary of the magnetic moments of the nucleons and light nuclei calculated with LQCD at 
SU(3) symmetric quark masses corresponding to a pion mass of m n ~ 806 MeV. The results are presented in 
units of natural nuclear magnetons. The red dashed lines correspond to the experimental magnetic moments. 


neutron star to lower its energy by spontaneously generating a large magnetic field. Second, a large 
isovector component to the nucleon magnetic polarizability is found. The proton polarizability is 
found to be considerably larger than that in nature while the neutron polarizability is consistent 
with the phenomenological value, but much more precise. Third, analysis of the j z = I z = 0 np 
system leads to a precise extraction of the coefficient, L\, of the short-distance two-body mag¬ 
netic current operator connecting the 3 S\ and \Sq states in the context of EFT(/). This operator 
provides an important contribution to the np —> dy capture cross-section near threshold, which is 
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FIG. 23: A summary of the magnetic polarizabilities of the nucleons and light nuclei calculated with 
LQCD at a pion mass of ~ 806 MeV. The upper panel presents the dimensionless quantity /3 = 
A— Mjv)/3/e 2 obtained from the fits with the inner shaded region representing the total uncertainty 
arising from statistical and fitting systematic uncertainties. The outer shaded region assesses additional 
systematic uncertainties from discretization effects and FV effects, combined in quadrature and applied 
multiplicatively. The lower panel presents the polarizabilities in physical units; in this case, the outer 
shaded region also includes the effect of the scale setting uncertainty. 
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a critical input for calculations of the production of elements in big-bang nucleosynthesis and in 
other environments as is discussed further in Ref [29) . 

These calculations are the first of their kind and are the initial steps in a comprehensive program 
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to determine the electromagnetic properties of the light nuclei as well as the response of nuclei to 
electroweak currents. The next steps will include calculations of axial matrix elements in the various 
light nuclei, as these are of significant phenomenological interest in neutrino-nucleus scattering 
experiments. Further, this points the way to calculating matrix elements of interactions required 
in laboratory searches for dark matter and other potential beyond the Standard Model scenarios 
which involve nuclear matrix elements of a variety of currents. Calculations at smaller lattice 
spacings and in other volumes, as well as for lighter quark masses where direct connection to 
experiment can be made are important to this program. Finally it is important to include the 
presently omitted couplings of sea quarks to the background fields. Calculations addressing these 
goals are planned for the future. 
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Appendix A: Charged-Particle Correlation Functions: Source Location and Gauge Origin 

1. General discussion 

In addition to a uniform magnetic field, the Abelian gauge links in Eq. lead to two further 
gauge-invariant quantities that are finite volume artifacts. These quantities are the Wilson loops 
Wi(x 2 ) and W 2 {xi) appearing in Eq. which express the non-vanishing holonomies of the 
background gauge field. A major consequence of these non-vanishing holonomies is the breaking of 
discrete translational invariance down to a smaller subgroup, which is referred to as the magnetic 
translation group, see Ref. [SO]. The size of this subgroup depends on the magnetic-flux quantum, 
h, as 


Wi 




1Ui(x 2 ), 


and 


W 2 




W 2 {x i) , 


(Al) 
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for j = 0,1,..., 3\q q \n — 1. Consequently lattice translational invariance for down-type quarks 
coupled to u\p\x), for example, is reduced to Zj x x Ll. This is to be contrasted with the 

a 

infinite-volume case, where full translational invariance in a uniform magnetic field is maintained 
due to gauge invariance. On a torus, gauge invariance is more restrictive due to the additional 
specification of Wilson loops, and translational invariance is consequently reduced. 

For charged-particle correlation functions, the reduced translational invariance can lead to subtle 
effects. For example, consider electromagnetic gauge links that are unity at the origin, ujp\ 0) = 1, 
such as those in Eq. Q, and the hadronic source to be at the spatial position X{. The electromag¬ 
netic gauge links can be altered so that they become unity at the source location, however, this 
cannot be achieved by a gauge transformation because the Wilson loops would be modified by 13 

Wi(z 2 ) -> and W 2 ( Xl ) ->■ W 2 (x 1 )W^(x iil ) . (A2) 

Similarly, translational invariance cannot be used to relocate the hadronic source to the origin 
without altering the correlation function. Consequently charged-particle correlation functions de¬ 
pend on the the origin of the gauge potential, and the location of the source. Even when Xi is 
related to the origin of the gauge potential by a discrete magnetic translation, which corresponds 
to the special case where the required translation is equivalent to a gauge transformation, the 
charged-particle correlation function will not be identical due to gauge dependence. 14 As the de¬ 
gree to which lattice translational invariance is reduced depends on the strength of the magnetic 
field, steps are required to ameliorate this situation. 

One way to deal with the problem is to fix the background held entirely including the holonomies. 
This approach can be implemented by randomly choosing a source location before including the 
ujp\x) links. The hadron source location must still be chosen relative to the gauge potential, and a 
convenient choice is to make them coincident. This method corresponds to implementing periodic¬ 
boundary conditions (BCs) on the quarks, 15 and was employed in |43] to investigate magnetic 
properties of the nucleon. The choice of a coincident location, however, is not required; and, other 
choices for the relative separation that are not related by a magnetic translation could be averaged 
over to mitigate FV effects. 

An alternative approach to reduce FV effects consists of varying the holonomies. One way to 
achieve this consists of introducing a constant shift in the gauge potential, which corresponds to 
implementing twisted-BCs on the quarks. 16 These are flavor dependent BCs due to the difference 
in quark electric charges. Ultimately to remove the arbitrariness of this choice, all non-equivalent 
shifts should be averaged over. The resulting twist average removes the FV effect associated with 
translational invariance, and related proposals have recently been suggested more generally to 
reduce FV effects in other lattice QCD computations [8T[ f82| . 17 


13 The electromagnetic gauge links are accordingly modified in the form —X [(a; — Xi) mod FJ. The 

new links U^ [(x — Xi) mod are related to Uj^\x — Xi) through a gauge transformation, however, such 
a transformation alters the functional dependence of charged-particle correlation functions. Because the Wilson 
loops are gauge invariant, there is no difference, for example, between Wi \{x2 — 2) mod and W\(x2—Xi t 2) = 

Wi(x2)W{(xi : 2), which appears above. 

14 The definition of the charged particle two-point function could be altered to include an electromagnetic gauge 
link between the source and sink locations. The resulting correlation functions would be gauge invariant; but, the 
cost is the introduction of a path for the gauge link. Because magnetic flux threads closed loops that are oriented 
transverse to the magnetic field, the resulting correlation function is then path dependent. 

15 Technically the quarks are only periodic (or twisted) up to a gauge transformation. This is often referred to as a 
magnetic periodic BC. 

16 Shifting the source and shifting the gauge potential are equivalent up to a gauge transformation, however, this 
leads to different expectations for charged particle correlation functions, as shown below. 

17 In the analogous case of time-varying gauge potentials which lead to electric fields, a variant of this procedure 
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The present calculations were performed with the following approach to handle the lack of 
lattice translational invariance. The source locations were varied relative to the origin of the 
gauge potential, rather than varying the gauge holonomies. Varying the source locations allows 
for the restoration of lattice translational invariance in two equivalent ways: averaging over all 
sources on a given configuration, and then performing the ensemble average; or performing the 
ensemble average with a fixed source location, and then averaging over all locations. In this 
work, a hybrid approach was chosen because of limited computational resources. Each QCD gauge 
configuration was post-multiplied by the Abelian gauge links uj^\x), and then a random offset 
was introduced. The random offset was the same for each value of the magnetic flux quantum, n, 
in order to maximize correlations between differing field strengths. On each configuration, quark 
propagators were calculated using 48 symmetrically distributed source locations. In addition to 
improving statistics, the source averaging partially restores lattice translational invariance on each 
configuration. Translational invariance is further improved in the ensemble average, because a 
random offset is chosen for each configuration. 

While quark propagators are subject to finite-size effects owing to the reduction in translational 
invariance, hadronic correlation functions for neutral particles are less susceptible. Using the neu¬ 
tron as an example, the non-vanishing holonomies of the gauge field show up only in exponentially 
small FV corrections to the neutron energy [84] . These corrections arise primarily from charged 
pion fluctuations that wind around the torus. For charged particles, however, the reduction in 
translational invariance has a direct effect on the coordinate dependence of their correlation func¬ 
tions. In turn, the overlap of a given hadronic interpolating field onto Landau levels depends on 
the hadron source location and gauge holonomy in a field-dependent way. To make this discussion 
more concrete, the simplified case of a point-like charged particle subject to the same method 
utilized in the present LQCD calculation is explored. 


2. Expectations for a Point-Like Charged Particle 

To illustrate the dependence of charged particle correlation functions on the source location 
and gauge holonomy, a point-like particle propagating in an external magnetic field on a torus is 
considered. In the point-like approximation, the spin-projected hadron propagator can be derived 
following the arguments in |56l [84] . Various ways to reduce finite-size effects by shifting the gauge 
potential and/or shifting the source location are considered. 

In practice, the point-like approximation is valid only when the typical hadronic size cannot be 
resolved in a given Landau level. For example, higher-lying Landau levels have a larger (rms) radius, 

oc yj ( til + \) /\Qh.eB\] hence, the details of the hadronic state will be less relevant compared with 
lower Landau levels. By contrast, the lowest Landau level is the most sensitive to hadron structure, 
and the most likely to be dynamically altered away from the point-like result. With magnetic fields 
that are not arbitrarily weak, more dependence on hadronic structure can be expected, and the 
point-like particle case may thus only provide a guide. Further study is needed to address the 
dynamical Landau levels of bound states, and to design better interpolating operators for hadrons 
in magnetic fields. 

Employing a uniform shift of the gauge potential transverse to the magnetic field direction, the 


was carried out for the neutron [ 83 ] • I 11 that study, results at second order in the gauge potential, were 
directly isolated by perturbatively expanding external field correlation functions, and enforcing temporal Dirichlet 
BCs. While twist averaging does not eliminate the FV effect in that case, it was shown that the neutron electric 
polarizability can nonetheless be separated from finite-size effects by efficacious shifts of the gauge potential. 
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gauge links are modified to 


u[ Q \x ) = exp[—x < 1 . ^° r Xl ^ L a 

I exp [—iqgn^-jj 2 - J for x\ = L — a 


U^\x) = exp[— iq q 02] exp 
U^\x) = 1 . 


iq q n- 


Q-Kax\ 

~1X~ 


The propagator for a structureless charged particle then has the form 

C hdz (x,x i; 0,B) = J2e- iQhe ^[wk^W 2 C^]{x + uL, Xi -B) , 


(A3) 


(A4) 


where x and Xi are the spacetime locations of the sink and source relative to the origin, respectively, 
and v is a triplet of integers . Images transverse to the magnetic field direction, = (zq,^), 
pick up phases arising from the constant shift of the gauge potential. Notice the Euclidean time 
direction is treated as infinite in extent. The propagator C^°^(x,xf, B) is that for a charged 
particle in infinite volume. Ignoring discretization effects, the continuum form of the infinite 
volume propagator is employed, which in coordinate space is |i8l ] 


C^°}(x,Xi',B) = W^(x,xi)^°](x - xf,B), 


h,jz 


(A5) 


and consists of two parts. There is a spacetime translationally invariant part 


(A 00 ) 

^h,j z 



ds QheBs 
(47 ts ) 2 smh.(QheBs) 



Q h eB 

Atanh(QheBs ) 


[zf + z 2 2 ] 


j- s [4 + z i\}, 

(A6) 


which contains the hadron energy Ehj z appearing in Eq. ([ 7 ]), however, the tilde denotes that 
it excludes the contribution to the energy from the n^-th Landau level. Instead, contributions 
from all Landau levels are contained in this propagator |85| . The remaining part (x, Xi) is not 
translationally invariant, and accordingly depends on the gauge. It can be written as a Wilson line 
evaluated on the straight-line path between Xi and x 


2U (x, Xi) = exp 


iQ h e / dz^AJz) 


= exp 


~-QheB(x 1 - x ii i)(x 2 + x i)2 ) 


(A7) 


In writing the FV propagator in Eq. (A4), the Wilson loops have been implicitly modified to 
include the hadron’s electric charge, Q^e, instead of the quark charge, q q e , and the charged hadron 
propagator C\j z (x,Xi;9,B) satisfies the following BCs in the directions transverse to the magnetic 
field: 


c h,j z {x + Lxi,Xi]6, B) = e lQhei C h .j z (x,Xi\6,B), 
Ch,j z (x + Lx2,Xi\0,B) = e lQh62 W2(xi)C h j z (x,Xi;0,B) 


(AS) 


The first is a twisted-BC, while the second emerges as a magnetic-twisted-BC. In the actual com¬ 
putation, the quark propagators are calculated with periodic BCs, however, the implementation of 
the external field with links that are periodic up to a gauge transformation, and on configurations 


with a gauge shift leads to the BCs in Eq. (A8). Because of the gauge shift, the gauge potential 


no longer vanishes at the origin, e.g., the links U^\x) are unity when x\ = O2L 2 /Quart. Finally, 


notice the asymmetric appearance of the holonomy W{x 1 ) in Eq. (A4). The Wilson loop VE(x 2 ) 
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does not appear explicitly in the sum over the winding number v \, and the charged-particle prop¬ 
agator is twisted in the aq-direction rather than magnetic twisted, see Eq. (A8). The effect of this 
Wilson loop, however, is contained implicitly in the sum over winding number iq because of the 
^-coordinate dependence of the Wilson line 2IJ(.T,a:j). This asymmetry in the charged-particle 
propagator and the BCs that result from it are directly related to the asymmetric choice of gauge. 

Given the form of the propagator in Eq. (A4), a natural question to ask is whether shifting both 
the gauge potential and the source location is superfluous. To answer this question, one can express 
the propagator in terms of the source-sink separation, Ax = x — x*, and attempt to absorb the 
remaining dependence on the source location into a redefinition of the twist angles 6. Due to the 
breaking of translational invariance, this is not possible. By virtue of the Wilson line 2B(x, Xi ), the 
correlation function retains explicit dependence on Xi 2 , which is measured relative to the origin. 
The origin has no significance for gauge-invariant quantities, however, in terms of the gauge links, 
the origin can readily be discerned, see Eq. (A3). As a consequence, gauge variant quantities, such 
as the charged-particle propagator, can depend on positions relative to the origin. 

Four scenarios are considered, denoted by T: i) periodic BCs with coincident origins of the 
gauge potential and source for the correlation functions (T = 0), ii) shifting the gauge potential 
(T = 6 ), iii) shifting the source location (T = X), and iv) varying the shift in the gauge potential 
and shifting the source location (T = OX ). 


a. Periodic BCs with coincident origins (T = 0) 


Choosing the origin of the gauge potential to coincide with that of the source and not including 
a uniform gauge field corresponds to specifying Xj = 0, and 6 = 0. The latter leads to periodic 
quarks (up to a gauge transformation). This method was chosen in [43] . Additionally in that 
study, as in this work too, the spatial sink location is summed over, which projects the correlator 
onto vanishing three-momentum. Three-momentum states do not have definite energy eigenvalues, 
however, and one expects correlation functions to contain multiple low-lying Landau levels. For a 
point-like particle on a continuous torus, consider the spatially-integrated correlation function, 

C 0 (t) = [ dx C h ,j z (x,0;0,B) . (A9) 

Jo 

Carrying out the three-momentum projection gives, 


C 0 (t) = 


1 r 

2 Jo 


ds 


_l s ®2 _ 1 t 2 

g 2 h,jz 4 s 


7TS 


pL 00 

/ dx 2 V] (x 2 , s\u 2 L, 0) , (A10) 

Jo 


v 2 =—oo 


which has been written in terms of the quantum-mechanical propagator for the simple harmonic 
oscillator 


(x',t'\x,t) = 6(t' — f)i 


Q h eB 


x exp < — 


2n sinh[Q/,ei3(t / — t)] 

QheB 


2sinh[Qf l eB(t' — t)] 


[(a/ 2 + x 2 ) cosh [QheB(t r — f)] — 2x'x\ > , (All) 


where t and t' are Euclidean times. In terms of Landau levels, the correlation function can be 
written as 


c 0 (t) = 4 ? 

n L =0 


e ^h,j z t 

^Bh,j z 


(A12) 
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where the energies Ehj, include the Landau energy, and are given in Eq. ([7]). The dimensionless 
coefficients Zn'l are given by 


l‘L °° 

Z n} = / dx 2 1p*n L {x 2) ^ 1 p nL (v 2 L) , (A13) 

J 0 V2=-oo 

and are written in terms of the coordinate wave functions, which have the standard form in terms 
of Hermite polynomials 

nL (x) = — - 1 H nL U\Q^B\x) e -HQ^ 2 . (A14) 

f nLnL -^B\ 

Notice that contributions from odd-parity Landau levels are absent due to the sum over winding 
number. The coefficients are not positive; spectral positivity is not maintained due to the lack 
of translational invariance. For quantized values of the magnetic field, these coefficients depend 
on the flux quantum h, but not on the size L. This dependence on n, however, is exponentially 
suppressed. 


b. Shifting the gauge potential (T = 8) 

While shifting the gauge potential has not been pursued as a means to reduce FV effects, it is 
instructive to discuss briefly point-like expectations for this method. Averaging over all possible 
shifts of the gauge potential is equivalent to averaging over quarks with randomly twisted boundary 


conditions. As a result, the twist-averaged propagator in Eq. (A4) receives contributions only from 


zero winding numbers = 0. While this is a desirable feature, there is no further simplification of 
the charged particle correlation function. The twists can be utilized, however, to form the infinite 
volume propagator via constructing the Fourier transformation in blocks |82| . The lack of lattice 
translational invariance means that the magnetic periodic images, rather than periodic images, 
must be summed over. In effect, this provides access to xf, B), for x± over the whole 

transverse plane. The resulting zero-momentum correlation function has the form 


Cg(t) = 


1 r c 

2 Jo 


ds 


— ^sE? . —f-t 2 

£ 2 h,j z 4 s 


y/2ns 


dx 2 (x 2 , s\x i} 2, o) 


(A15) 


In terms of Landau levels, the expected behavior is thus 


Ce(t) = J2 zi n] 


n L =0 


2 E k 


(A16) 


Jz 


where the coefficients are given by 

/ oo 

dx 2 lf* nL {x2)^n L {Xi,2) , (A17) 

-oo 

in which there is no remaining dependence on L. When the source is located at the origin, there 
is no dependence on the magnetic flux quantum, n. Notice that there are no contributions from 
odd-parity Landau levels, and further that spectral positivity does not emerge. 
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c. Shifting the source location (F = X) 


As described previously, source-to-all propagators are calculated in the present LQCD calcu¬ 
lation, with sources randomly located relative to the gauge origin. A constant shift of the gauge 
potential is not implemented, and thus 0 = 0 , but with an approximate average over X{. This 
leads to the following expression for the correlation function of a point-like charged particle, 


Cx(t ) 


1 

2 


ds - 


-lsE 2 h . - 

- ° h,Jz 



(X 2 ,S \x ij2 ,0) 


(A18) 


Performing the integration over source location restricted the winding number expansion to 
the sector with v -2 = 0. Using the spectral decomposition for the quantum-mechanical harmonic- 
oscillator propagator gives, 


c x (t) = J2 4 ? 

n L =0 



(A19) 


where Ehj z includes the energy of the n^-th Landau level. The corresponding spectral weights are 


Z™ = - 
L 


dx lf nL (x) 


(A20) 


and give the probability of finding the charged particle in the n^-tli Landau level. Positivity of these 
weights arises due to the symmetric treatment of the source and sink locations. Landau levels of 
both parities contribute to the correlation function. When evaluated for quantized magnetic fields, 
the weights only depend on the flux quantum n. Ratios of coefficients, however, depend on h 
through exponentially small terms. 


d. Varying the gauge shift and source location (r = OX) 

Despite computational requirements, it remains worthwhile to consider averaging over the shift 
of the gauge potential and the source location. The former can be used to construct magnetic- 
periodic images and build the infinite-volume propagator in blocks. The result of this procedure 
leads to 


Cexit) = 2 

n L =0 



(A21) 


where Ef l j z includes the energy of the n^-th Landau level, and the corresponding spectral weights 
are 


Z& x ) 

n L 


1 

L 


dx 'lp rlL (x) 


(A22) 


This procedure eliminates all L dependence, and n dependence is completely absent in ratios of 
coefficients. Furthermore, the procedure excludes contributions from odd-parity Landau levels, 
and maintains spectral positivity. 
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3. Results for a point-like particle 


Having determined the charged-particle correlation functions for four different methods of deal¬ 
ing with the lack of translational invariance, the relative contributions of the lowest-lying Landau 
levels are compared in Fig. 24 In plotting the coefficient ratios, the n^-dependence of the hadron 
energies Ehj z appearing in the correlation functions has been ignored. These energy denominators 
will lead to smaller coefficients, but only for higher-lying Landau levels. 



n L 

FIG. 24: Contributions to the correlation functions of a point-like charged particle from the first few 
Landau levels compared with the contribution from the lowest Landau level. Four such coefficient ratios 
are considered, corresponding to different ways of dealing with the lack of lattice translational invariance, 
as described in the text. These are labeled by T = 0, 6, X , and OX , which correspond to: coincident origin, 
twist-averaged, source-averaged, and twist and source-averaged, respectively. Values have been slightly 
displaced in to allow different ratios to be discernible. For quantized magnetic fields, all ratios are 
independent of the lattice size L , and ratios are either independent of the flux quantum n, or depend on it 
only through exponentially small terms. 
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